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C ~ 3 ■ We review recent progress in understanding QCD at high energies and the role 

' played in it by the large JV C limit. We discuss unitarization of total hadronic cross 

^ _ sections and saturation of structure functions at high energies. 

1 The BFKL Equation 

in : 

One of the most important problems in QCD is the problem of understand- 
ing its high energy asymptotics. High energy QCD has direct relevance to 
the scattering data produced at the accelerators worldwide allowing for ex- 
perimental verification of theoretical attempts to understand QCD dynamics. 
C*") ■ In addition to that gluonic fields in this extreme regime of high energies be- 

come very strong making the underlying description of QCD processes non- 
perturbative, while leaving the strong coupling constant a s small. Therefore 
\ understanding the gluon dynamics and correct degrees of freedom of QCD 

■ at high energies is possible and may shed some light on the non-perturbative 
dynamics of the theory. 

• Let us consider a scattering of two heavy quarkonia (quark-antiquark 

1 \ states), or, equivalently, scattering of two virtual photons with large virtu- 

alities Q\ and Q\, each of which splits into a quark-antiquark pair. Large 
(-h \ virtualities or large masses of the quarkonia insure that the strong coupling 

L" ■ constant is small, a s {Q\) ~ ct s (Q\) <C 1, allowing us to use perturbative 

. £^ \ expansion. 

■ At the lowest order in a s the simplest model of high energy interaction 
between the two quarkonia is given by a two gluon exchange in the forward 
amplitude. The corresponding total cross section of onium-onium scattering 
is independent of energy at very high energies and could be written as 

(1) 

where s is the center of mass energy of the system. This model of high 
energy interactions is known as Low-Nussinov pomeron.EI It demonstrates 
that exchanging two vector particles (gluons) between the onia gives a non- 
vanishing cross section at high energy, while using scalars or fermions (with 
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Figure 1. Ladder diagrams summed by the BFKL evolution equation. 



non-derivative couplings) would give a cross section that decreases with en- 
ergy. 

However, at the high energies achievable by the modern days accelerators 
the hadronic cross sections are not constant. They tend to increase as powers 
of center of mass energy of the system 



where A is a number called pomeron intercept. The rise of partonic structure 
functions with energy could be described by the DGLAP evolution equationB 
Nevertheless the DGLAP equation does not yield us with a power of energy 
growth of the cross sections. The only way to obtain a cross section in pertur- 
bative QCD that grows as a power of energy is given by the BFKL evolution 
equations n 

The BFKL equation a resums ladder diagrams in the forward onium- 
onium scattering amplitude of the type shown in Fig. [l| The light cone 
momenta components of the s-channel gluons are ordered (Regge kinematics) 
such that 



(2) 



p+ > k 1+ > k 2+ ...» k N+ > p' + 



(3) 



and 



p- < k X - < k 2 - ■ ■ ■ < k N _ < p'_ 



(4) 
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while there is no restrictions on their transverse momenta which remain of 
the same order throughout the ladder 

Q 1 ~ k\x ~ k 2 ± ~ . . . ~ k N± ~ Q 2 . (5) 

Each s-channel gluon in the ladder brings in an extra power of the coupling 
constant a, « 1. The integration over the longitudinal phase space of each 
gluon gives a term proportional to In s 3> 1 . Combining the two factors we 
see that each extra s-channel gluon introduces a factor of a s Ins ~ 1. Thus, 
by resumming ladders of the type shown in Fig. |l| the BFKL equation resums 
all powers of the effective parameter a s Ins. 

By drawing a ladder in Fig. [l] we are oversimplifying the problem. The 
triple gluon vertices in the ladder are not the usual QCD triple gluon vertices. 
These are the so-called effective Lipatov vertices and have to be calculated by 
analyzing several gluon emission diagramsH, though in certain gauges they do 
correspond to the usual triple gluon vertices. The gluon propagators in the 
t-channel are also not just the usual gluon propagators. They include virtual 
(reggeization) corrections to the ladder which modify the gluon propagators. 
The resulting reggeized gluons are denoted by thick lines in Fig. each line 
representing a summation of the whole ladder of Fig. [I] in the color octet 
channel. 

The BFKL equatinn-jCould be written for onium-onium scattering ampli- 
tude f(Qi,Q 2 ,Y) asW&B 



df(l,L',Y) 
dY 



= J d 2 kK BFKL (l_,k)f(k,!!,Y) (6) 



where Y — ln(s/Qi<52) is the rapidity variable and Kbfkl is the kernel of 
the equation. Defining the Laplace-transformed amplitude as 

/>oo 

fM) = / dYe-» Y f(l,l',Y) (7) 



Eq. (H) could be rewritten as 

aN r f d 2 k 



ufM) = s 2 (i_-i!) + 



fM')- l2fJl - ] 



k 2 + {k- If 



(8) 



The solution of Eq. (||) in the saddle point approximation at large Y yields 
the scattering amplitude 



f(Qi,Q2,Y) = e (^-!) y exp 

2nQ 1 Q 2 V^DY 



ln 2 Qi/Q 2 



4DY 



(9) 
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with D = ^iV c C(3) and a P - 1 = ^^ln2. At very high energy the 
solution given by Eq. (P) grows as a power of energy since 



,(ap— l)Y ap — l 



(10) 



Therefore the BFKL equation seems to be a natural candidate for the phe- 
nomenological pomeron as observed in high energy scattering experiments. It 
is beyond the scope of this talk to address the issue of how well the BFKL 
equation describes the data. We just note that the numerical value of the 
obtained pomeron intercept ap — 1 = ias J sr " In 2 ps 2.6 a s « 0.8 is higher 
than the experimental value of 0.2 — 0.3 observed in deep inelastic scattering 
experiments. The inclusion of next-to-leading order correction to.JBFKL ker- 
nel makes the intercept negative for all reasonable values of a s B0l3, though 
this problem is likely tq-.be due to collinear singularities and should be cured 
at higher orders in a s □ resulting in a slightly lower but still positive value 
of the intercept which is closer to experimental data. The BFKL equation 
as described here is applicable only to the scattering of perturbative (small) 
objects like heavy quarkonia and can not be applied to the total cross section 
in, say, proton-proton collisions, where (in the single ladder approximation) 
there is no hard scale which would justify the small coupling approach. 

However, the BFKL equation poses some important questions even in the 
case of heavy onium-onium scattering and at the leading order in a s in the 
kernel. 

(i) The power of energy growth of the total cross section (Eq. (|l(])) vio- 
lates Froissart unitarity bound, which states that the growth of the total cross 
sections with energy at asymptotically high energies is bounded by 113 

const 2 

a < 77- In s (11) 

ml 

with m v the pion's mass. XFor a good pedagogical derivation of Froissart 
bound I refer the readers to llil.) That means some new effects should modify 
the BFKL equation at very large s to make the resulting amplitude unitary. 

(ii) The solution of Eq. (g) includes a diffusion term. To see this let 
us consider a half of the ladder, from one of the onia to some intermediate 
gluon in the middle of the ladder carrying momentum k±_ and rapidity Y/2. 
Applying Eq. <M) to that half-ladder we see that it includes a term 



exp 



In 2 k ± /Q 
2DY 



(12) 



This term is responsible for diffusion of the transverse momenta from the 
initial perturbative scale Q both into infrared and ultraviolet. The width of 
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the diffusion grows with rapidity Y. Thus no matter how large the starting 
scale Q is at certain very high energy the momentum of some gluons in the 
middle of the ladder would become of the order of A^qcd invalidating further 
application of BFKL evolution£3 Again this hints that the BFKL equation 
should be modified at higher energies to avoid this problem. 

2 Multiple Reggeon Exchanges 

The BFKL ladder in Fig. g can be viewed as two reggeized gluons (reggeons) 
propagating in the t-channel. It has been conjectured that unitarity of the 
total cross section could be restored if one resums all diagrams with multiple 
reggeon exchanges E3. Instead of summing the ladder in Fig. [I] one has to 
analyze all the diagrams of the type shown in Fig. ^. There one can exchange 
any number of reggeons in the t-channel. Each pair of reggeons may interact 
via s-channel gluon exchanges leading to BFKL-type contributions. To com- 
plete the picture one should also include the diagrams where the number of 
reggeons varies in the i-channel, that is the diagrams containing vertices of n 
reggeons going into m reggeons for any n and m. This is beyond the scope of 
this talk. We are going to concentrate on the diagrams where the number of 
reggeons is constant. 

The n-teggeon exchange diagrams were first resummed by Matinyan and 
Sedarkyan llil for the scalar </> 3 theory. For n-reggeon exchange amplitude 
one can write down an equation similar to Eq. M). However there are many 



Figure 2. Multiple reggeon exchange diagram. 
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more contributions to the kernel than in the 2-reggeon case considered above. 
Since at each step in rapidity any two out of n reggeons can interact one 
expects n(n — l)/2 contributions to the kernel of the integral equation. Each 
contribution comes in with the same integral kernel acting on a different pair 
of reggeons. Assuming for simplicity that all these kernels are approximately 
the same we can write the following equation for n-reggeon amplitude T n at 
large n 



dY 



n 2 



K®F n . (13) 



The solution of Eq. (13) scales as 

■p (Y'j g(-2n+COnstn 2 )Y ^ s -2n+COnstn 2 

where the term — 2n in the power results from the fact that we are dealing 
with n scalar particles in 3 theory. The total cross section would be given 
by the sum of n-reggeon amplitudes 

E-^ ~ E c « s ~ 2 " +const " 2 > ( 15 ) 

n n 

which is not likely to converge! (c„'s are some coefficients which may depend 
on s at most logarithmically.) This result indicates that multiple reggeon 
exchanges in <fi 3 theory not only fail to unitarize the total cross section, but 
rather make the series highly divergent. It is believed that the problem is 
associated to the fact that 3 theory has no stable vacuum and is typical only 
to 4> 3 theory. Nevertheless one might and should get worried about similar 
phenomenon happening in QCD. 

To demonstrate that the catastrophe which happened in <f) 3 theory is 
actually avoided in QCD one has to invoke the 't Hooft large N c limit t3 which 
is the central topic of this conference. In the large N c limit the diagrams in 
Fig. ^ drastically simplify: only the planar graphs survive. Thus the topology 
of the interaction takes form of a cylinder with the reggeons forming the walls 
of the cylinder and only the interactions of the nearest neighbor reggeons are 
allowed, as depicted in Fig. ^. There the ith reggeon may interact only with 
the i— 1st and i + 1st reggeons, and the 1st reggeon interacts with the nth and 
the 2nd reggeons. The two reggeon case (the pomeron of Fig. |l|) corresponds 
to the simplest cylinder topology case with reggeons and s-channel gluons 
"living" on the walls of a cylinder. If we add one more reggeon we would 
obtain the three reggeon case (odderon) which also has a cylindrical topology 
as shown in Fig. |^. 

With only the nearest neighbor interactions between the reggeons allowed 
the kernel of the integral equation resumming the diagrams shown in Fig. [| 
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Figure 3. An example of cylinder topology acquired by multiple reggeon exchange diagrams 
in the large N c limit. The ellipses are drawn to guide the eye. 



receives only n contributions, which is very different from (f> 3 theory. Making 
the same assumption as for <f> 3 theory we see that Eq. @ for large- N c QCD 
becomes 

QT 

-^^nK®T n (16) 

leading to 

T n {Y) ~ e constnY _ s constn (17) 
Now the total cross section similar to Eq. (JT^) is 

a „ J2?n - ^d„ S const ", (18) 

n n 

with d n 's some coefficients with weak dependence on s. The series in Eq. (HI 
could be convergent and therefore large- N c QCD avoids the catastrophe of < 
theory. One expects that the problem would not reappear for the real life case 
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of three colors. We have to advise the reader to understand Eq. ([L7|) only as 
an upper bound on the n reggeon amplitudes but not as an exact answer. 

A careful QCD analysis of the reggeon cylinder graphs of the type shown 
in Fig. U demonstrates that a reggeon system is equivalent to a one dimensional 
XXX Heisenberg spin chain for spin S — and is thus an exactly solvable 
model Ell. That facilitates the analysis of the n-reggeaa diagrams allowing to 
quantitatively understand their high energy behavior.113 

An important question for the total cross section is whether the leading- N c 
2n-reggeon diagram gives an amplitude which grows with energy faster than 
the subleading in N c n-pomeron diagram. Naturally the n-pomemn diagram 
should scale as s "( a J'- 1 ) . I n the recent paper by Korchemsky et alia it has been 
shown that the n reggeon amplitude seem to scale as T n ~ s const/n^ Thus 
at very high energies multiple reggeon exchange diagrams grow much slower 
than multiple pomeron exchange diagrams and the former could be neglected 
compared to the latter. We refer the interested readers to the contribution of 
G. Korchemsky in these proceedings 113. 



3 Multiple Pomeron Exchanges 

3.1 General Physical Picture 

As could be seen from the discussion in the previous section an exchange of a 
single BFKL pomeron brings in a parametrical factor of a 2 e ( a p~ 1 ) Y ; where 
e (ap-i)Y comes f rom the amplitude in Eq. (^) and the factor of a 2 comes 
from connecting the reggeons in the ladder to the onia. The BFKL equation 
resums powers of [ap — 1)Y ~ a s Y and therefore gives a large contribution 
when 

Ylo ~ — ■ (19) 

a s 

What happens as we increase the energy /rapidity? It is possible that mul- 
tiple pomeron exchanges would start playing a bigger role in the scatter- 
ing. To estimate when this happens we have to equate single (a 2 e^ ap ^ 1 " >Y ) 
and double (a 4 s e 2( - ap -^ Y ) pomeron contributions, ^effectively requiring that 
a 2 e (a P -i)Yu ^ i was ar g ue d by Mueller in 113 multiple pomeron ex- 
changes become important at the values of rapidity of the order of 

Yu ^ln^. (20) 

ap — 1 aj. 

After completion of the calculation of the next-to-leading order correc- 
tions to the kernel of the BFKL equation (NLO BFKL) by Fadin and Lipatov 
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Figure 4. An example of the pomeron "fan" diagram. Each ladder interacts with the proton 
or nucleus at rest independently. 

i, and, independently, by Camici and Ciafalonii, it was shown in I that due 
to the running coupling effects these corrections become important at the 
rapidities of the order of 

Ysjuo - 4fr ^ 

One can see that Ylo <C Yjj ^ Ynlo f° r parametrically small a s . That im- 
plies that the center of mass energy at which the multiple pomeron exchanges 
become important is much smaller, and, therefore, is easier to achieve, than 
the energy at which NLO corrections start playing an important role. That 
is multiple pomeron exchanges are probably more relevant than NLO BFKL 
to the description of current experiments.E3 Also multiple pomeron exchanges 
are much more likely to unitarize the total hadronic cross section. Here we 
are going to present a solution to the problem of resummation of multiple 
pomeron exchanges and show how the BFKL pomeron unitarizes followingrll. 

We will consider deep inelastic scattering (DIS) of a virtual photon on a 
hadron or nucleus and will resum all multiple pomeron exchanges contributing 
to the total DIS cross section in the leading longitudinal logarithmic approx- 
imation in the large N c limit. The first step in that directicp_in perturbative 
QCD was a conjecture by Gribov, Levin and Ryskin (GLR)c3 of an equation 
describing the fusion of two pomeron ladders into one, which was proven in the 
double logarithmic limit by Mueller and Qiu.E3 The resulting equation resums 
all pomeron "fan" diagrams (see Fig. ^) in the double logarithmic approxima- 
tion. An extensive work on resumming the multiple pomeron exchanges in the 
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Figure 5. Dipole evolution in the deep inelastic scattering process as pictured i nEI. The 
incoming virtual photon develops a system of color dipoles, each of which rescatters on the 
target hadron or nucleus (only two are shown). 



gluon distribution function in the leading ln(l/x) approximation (i.e. without 
taking the double logarithmic limit) both in the framework of effed.ipj.iicM 
theories and employing alternative approaches has been pursued in CJ'E30E3. 

Let us consider a deep inelastic scattering (DIS) of a virtual photon on a 
hadron or nucleus. An incoming photon splits into a quark-antiquark pair and 
then the qq pair rescatters on the target hadron or nucleus. In the rest frame 
of the hadron all QCD evolution should be included in the wave function of 
the incoming photon. That way the incoming photon develops a cascade of 
gluons, which then scatter on the hadron at rest. We want to calculate this 
gluon cascade in the leading longitudinal logarithmic (In 1 jx) approximation 
and in the large N c limit. This is exactly the type of cascade described by 
Mueller's dipole modeEa. In the large N c limit the incoming gluon develops 
a system of color dipoles and each of them independently rescatters on the 
hadron or nucleus. The forward amplitude of the process is shown in Fig. ^. 
The double lines in Fig. || correspond to gluons in large N c approximation 
being represented as consisting of a quark and an antiquark of different colors. 
The color dipoles are formed by a quark from one gluon and an antiquark from 
another gluon. In Fig. || each dipole, that was developed through the QCD 
evolution, later interacts with the hadron or nucleus by a series of Glauber- 
type multiple rescatterings on the sources of color charge. The assumption 
about the type of interaction of the dipoles with the hadron is not important 
for the evolution. It could also be just two gluon exchanges. The important 
assumption is that each dipole interacts with the target independently of the 
other dipoles, which is done in the spirit of the large N c limit and is valid for 
a model of a large dilute target nucleus. 

For the case of a dilute target independent dipole interactions with the 
target proton or nucleus are enhanced by some factors of the number of 
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sources of color charge in the proton or nucleus compared to the case when 
several dipoles interact with the same parton in the proton or the same nu- 
cleon in the nucleus .c3ti3 This allowed us to assume that the dipoles interact 
with the proton or nucleus independently. In case of a nucleus one actually 
has the atomic nurabei A as a parameter leading powers in which justify 
this approximation.cj'E2l The summation of these contributions corresponds 
to summation of the pomeron "fan" diagrams of Fig. [|. In general there is 
another class of diagrams, which we will call pomeron loop diagrams. In a 
pomeron loop diagram a pomeron first splits into two pomerons, just like in 
the fan diagram, but then the two pomerons merge back into one pomeron. 
The graphs containing pomerons not only splitting, but also merging together 
will be referred to as pomeron loop diagrams. In a dilute proton or nucleus 
these pomeron loop diagrams could be considered small, since they would be 
suppressed by powers of color charge density (A) . At very high energies this 
approximation becomes not very well justified even in the dilute target case. 
Contribution of each additional pomeron in a fan diagram on a hadron is 
parametrically of the order of a^A 1 ^ 3 e^ ap ~ 1 ^ Y which is large at the rapidities 
of the order of Yjj given by Eq. ((2^). The contribution of an extra pomeron 
in a pomeron loop diagram is of the order of a 2 s e lyaF ^ 1 ^ Y , which is not en- 
hanced by powers of A and is, therefore, suppressed in the dilute nucleus case. 
However, one can easily see that at extremely high energies, when a 2 e^ ap ~ 1 ' >Y 
becomes greater than one, pomeron loop diagrams become important, still be- 
ing smaller than the fan diagrams. Resummation of pomeron loop diagrams 
in a dipole wave function seems to be a very hard technical problem, possibly 
involving NLD BFKL or, even, next-to-next-to-leading order BFKL kernel 
calculations £il Nevertheless, it is the author's belief that while being crucial 
for onium-onium scattering the pomeron loop diagrams are not going to sig- 
nificantly change the high energy behavior of the DIS structure functions. 

3.2 Nonlinear Evolution Equation 

We start by considering a deep inelastic scattering process. The incoming 
virtual photon with a large q+ component of the momentum splits into a 
quark-antiquark pair which then interacts with the target at rest. We model 
the interaction by no more than two gluon exchanges between each source 
of color charge and the quark-antiquark pair. This is_done in the spirit of 
the quasi-classical approximation used previously in E3'LJ. The interactions 
are taken in the eikonal approximation. Then, as could be shown in general, 
i.e., including the leading logarithmic QCD evolution, the total cross-section, 
and, therefore, the F2 structure function of the nucleus can be rewritten as a 
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product of the square of the virtual photon's wave function and the propagator 
of the quark-antiquark pair through the proton. The expression reads 

F 2 (x,Q 2 )= - f j ^i^[$ T (x 01 ,z) + $ i (x 01 ,z)]d 2 &o^(xoi,b ,y), 
iTT z a E M J 2tt 

. (22) 

where the incoming photon with virtuality Q splits into a quark-antiquark 
pair with the transverse coordinates of the quark and antiquark being xo 
and xi correspondingly, such that xio = xi — xo. The coordinate of the 
center of the pair is given by bo = hfal + Xo). The square of the light cone 
wave function of qq fluctuations of a virtual photon is denoted by $t(xoi, z) 
and $i(xoi,vz) for transverse and longitudinal photons correspondingly, with 
z being the fraction of the photon's longitudinal momentum carried by the 
quark. (E > t(xoj- t1 2;) and $l(xoi, z) are known functions and could be found in 
the literatureO 

The quantity 7V(xoi,bo,F) has the meaning of the forward scattering 
amplitude of the quark-antiquark pair on a nucleus. At the lowest (classical) 
order not including the QCD evolution in rapidity it is given by 

-x^AxGfol/xox) 



iV(x O i,b o ,0) = -7(x i,b ) ~ < 1 - exp 



2N C S ± 



(23) 

7(x i , bo) is the propagator of the qq pair through theJiadron or nucleus. The 
propagator could be easily calculated, similarly to E3S, giving the Glauber 
multiple rescattering formula (p3|). Throughout the talk we assume for sim- 
plicity that the target nucleus is a cylinder, which appears as a circle of radius 
R in the transverse direction and has a constant length 2R along the longitu- 
dinal z direction. Therefore its transverse cross-sectional area is S± = ttR 2 . 
In formula ( p3|) a s is the strong coupling constant and xG(x,1/xq 1 ) is the 
gluon distribution of a nucleon in the nucleus, taken at the lowest order in a s , 
similarly to E3. 

Eq. ( p3| ) resums all Glauber type multiple rescatterings of a qq pair on 
a target nucleus. As was mentioned before, since each interaction of the 
pair with a source of color charge in the nucleus is restricted to the two gluon 
exchange, the formula ( |23| ) effectively sums up all the powers of the parameter 
a^A 1 / 3 . (Since xG ~ Cf there is no N c present in the expansion in the large- 
N c limit.) Therefore this interaction of course is subleading in N c as is true 
for any hadronic interactions.!!^ 

Since the target is at rest in order to include the QCD evolution of F 2 
structure function, we have to develop the soft gluon wave function of the 
incoming virtual photon. In the leading longitudinal logarithmic approxima- 
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Figure 6. A diagrammatic representation of the dipole evolution equation. 



tion (In 1/x) the evolution of the wave function is realized through successive 
emissions of small-a; gluons. The qq pair develops a cascade of gluons, which 
then scatter on the hadron. In order to describe the soft gluon cascade we 
will take the limit of a large number of colors, N c — > oo. Then, this leading 
logarithmic soft gluon wavefunction-will become equivalent to the dipole wave 
function, introduced by Mueller in l3. The physical picture becomes straight- 
forward. The qq pair develops a system of dipoles (dipole wave function), and 
each of the dipoles independently scatters on the target hadron or nucleus (see 
Fig. ||) . Since the target is assumed to have many sources of color charge in it 
we may approximate the interaction of a dipole (quark-antiquark pair) with 
it by 7(x, b) given by Eq. (p3|), with x and b being the dipole's transverse 
separation and impact parameter. That means that each of the dipoles inter- 
acts with several sources of color charge (Glauber rescattering) in the nucleus 
independent of other dipoles. 

To construct the_dipole wave function we will heavily rely on the tech- 
niques developed in c3. We have to resum a gluonic cascade in the large- N c 
limit. Denoting the cascade by a blob we can write down an equation which is 
illustrated in Fig. |(| As usual for the large N c limit the double line denotes a 
gluon. The original qq pair may have no evolution in it, i.e., the cascade could 
be empty, as illustrated by the first term on the right hand side of Fig. 0. In 
one step of the evolution in rapidity a gluon is emitted in a dipole, splitting 
the parent dipole into two dipoles. Each of the produced dipoles is nonlo- 
cal: it consists of a quark (antiquark) from the original quark-antiquark pair 
and an antiquark (quark) component of the emitted gluon. The consecutive 
gluon evolution can continue in either of the two dipoles. This is depicted in 
the second term on the right hand side of Fig. |[ This, procedure gives us a 
self-consistent equation resumming the gluon cascadeH3'Eil 

Substituting 1 — iV(xoi, bo, Y) instead of the blob symbolizing the gluon 



cascade and including virtual corrections one obtains the following equatio: 
iVfxoi.ba, Y) =-7(x i,b ) exp 
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ol s Cf 



dy exp 



d 2 X2 



x oi 

2 2 
a: 02 :Z: 12 



[2JV(xo2,bo + -xi2,y) 



7V(xo2,bo + ixi2,y)iV(xi2,bo - ^x 20 ,y)]. 



(24) 



Equation (J2J), together with Eq. (E22h, provide us with the leading logarith- 
mic evolution of the F 2 structure function of a hadron including all multiple 
pomeron exchanges in the large- N c limit. Note that while the evolution is 
leading in N c the interactions (7) are not, so that Eq. ( pi]) could be viewed 
as resumming powers of a 2 A 1 / 3 e °>sN c Y _ 

Similar evolution equation was obtained by Balitsky in c3 using an ef- 
fective lagrajagian approach to high energy QCD. The evolution equations 
obtained in Ei3 reduce to a closed sinadfi|-epuation only in the large- N c limit. 
Similar results have been obtained in rirl. 



3.3 Solution of the Nonlinear Evolution Equation 



To understand the qualitative features of the integral equation ( 24 ) it is con- 
venient to rewrite it in the differential form. Assuming that the typical size 
of the dipole wave function x± <C R, where R is the proton's radius, we can 
rewrite Eq. (p4j) as 



cW(x i,F) 2a s C F , , 2 

d x 2 



dY 



'01 



2 2 



27T(5 2 (xoi - x 02 )ln 



Ot s CF 



7T J X Q2 X 12 

Performing a Fourier transform 



d'x 2 iV( X0 2, F) JV(xia, Y). 



N(x 2,Y) 
(25) 



N(x±,Y) = x\ 



~2T 



e 4k - x N(k, Y)=x± dkk J (kx±) N(k, Y). (26) 



we obtain 

dN(k,Y) 2a s N c 



dY 



d 



d\nk 



N ( ki Y)-^^N 2 (k,Y), (27) 



where 



(28) 
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"3>- 



energy s 



Figure 7. Qualitative behavior of the total DIS cross section given by the solution of 
Eq. (g). 



is the eigenvalue of the BFKL kernel I with ip(\) = T'(X)/T(X). In Eq. 
the function x(A) is taken as a differential operator with A = — <9/<91nfc acting 
on N(k, Y). We also put N c /2 instead of Cp everywhere in the spirit of the 
large iV c approximation. The details of obtaining Eq. (|2^) from Eq. (p5| ) 
could be found in Ell. 

Now the physical meaning of the evolution given by Eq. (|24|) becomes 
manifest. The linear part of Eq. (|27| ) is just the familiar BFKL equation, 
giving a cross section growing as a power of energy. When energy is not too 
large and N is small the quadratic term could be neglected and the linear 
evolution is all that remains. However as the energy increases the evolution 
driven by BFKL would make N large, making the quadratic term in Eq. ( p7| ) 
important. The quadratic term would introduce damping effects, similar to 
GLR-MQ equationE§E3, slowing down the growth of the cross section and 
possibly unitarizing it. r - ir - 1 

Eq. ( p4| ) has^eHn-snlved both analytically by approximate methods E1E3 
and numericallyE§L3EJ'E3. The result for DIS cross sections is qualitatively 
illustrated in Fig. pi The cross section starts growing as-a-japwer of enerj 

BET 



and then saturates to a much slower (logarithmic) growth. BIB Thus Eq. 
resolves the problem of unitarization of BFKL pomeron (i) posed in Sect. 1. 



kovchegov: submitted to World Scientific on February 1, 2008 15 



Does Eq. (p3) answer the second question (ii) posed in Sect. 1? Let 
us first note that Eq. ( p^ ) produces a large momentum scale which grows 
with enemy- as predicted within the framework of McLerran-Venugopalan 
modeSEflEa. To see this let us write down the solution to Eq. (pTj) in the 
limit of moderate energies (where the equation is still mostly linear): 

~ A exp[(« P - 1)Y] / ZL__ *) , (29) 

where C is some overall normalization constant and A is the momentum scale 
inherent to the initial distribution given by Eq. (^3|) (in—facl-it is the classical 
saturation scale given by McLerran-Venugopalan model_rE!l_3) . At very high 
energy the amplitude N saturates tp_a constant: N — > 1 as Y — > oo. In 
momentum space this corresponds tot_J 

N sat (k,Y) =ln^. (30) 

Logarithm is a slowly varying function compared to exponent and to estimate 
at what momentum scale the transition to saturation takes place we can just 
equate the amplitude in Eq. ( ^9| ) to one. Omitting all slowly varying with Y 
prefactors we obtain an estimate of the the saturation scale 

Q a (Y) w A BX ^ ap ~ 1)y] . (31) 
^/Ua s N c C(3)Y 

One can see from Eq. (|3|) that as energy increases the saturation scale in- 
creases too. j_.j_.j_. 

Numerical simulation— l'—l't_ demonstrated a very interesting behavior of 
the solution of Eq. (pi|), which is drastically different from the solution of 
BFKL. Let us remind the readers that in Sect. 1 based on Eq. ( [12] ) we ob- 
served that with increasing energy the momentum of gluons given by the 
solution to BFKL equation diffuses around some initial external perturbative 
scale Q with the diffusion width increasing with Y. That was potentially dan- 
gerous since the diffusion cone would eventually get into the non-perturbative 
region. The solution of Eq. ( |2^ ) produces a very different distribution of 
transverse momenta of the gluons: the width of the diffusion cone seems to 

.tribution is given by the 
I (seeEq. @). Thus the 



be independent of energy while the center of ilu 
saturation scale Q S {Y) which increases with Y_§ 

distribution never gets into the dangerous infrared region. Therefore Eq 
appears to solve the problem (ii) posed in Sect. 1 of the sensitivity of the cd 
sections given by BFKL evolution to the non-perturbative infrared region.E. 
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4 "Phase Diagram" of High Energy QCD 

To summarize the results of the previous Section we present the "phase dia- 
gram" of high energy QCD depicted in Fig. ||. 

The majority of the HERA data shows that the proton's structure func- 
tions at large photon's virtuality Q 2 and not very^mall Bjorken x w Q 2 / s can 
be described by the DGLAP evolution equation □, which is a linear equation. 
It is convenient to present different properties of hadronic structure functions 
in terms of a "phase diagram" in Q 2 and lnl/x plane (see Fig. ^). The 
DGLAP physics corresponds to the lower right section of the plane, where Q 2 
is large and x is not too small. As one goes towards smaller x in the same 
region of Q 2 the hadronic structure functions rise. Most of the data in that 
region can be explained in terms of cither small- a; limit of DGLAP equation 
or, alternatively and more interestingly, by the BFKL equationtj, which is also 
a linear equation but could be responsible for evolving the system towards a 
higher gluonic density regime. However, the DIS cross sections can not rise 
forever as powers of center of mass energy s. This would violate unitaritytL9 
and pose other problems discussed in Sect. 1. That means that at some 
very small x the hadronic structure functions have to undergo a significant 
qualitative change of behavior, becoming a much slower varying functions of 
x. The slow down of the growth of hadronic structure functions is usually 
associated with non-linear effects in the quark and gluon dynamics, such as 
parton recombination, which eventually balances the parton splitting process 



In 1/x 



non-perturbative 
region 
(not much is known 
coupling is large) 
a ~ 1 



saturation region 

( can be understood 
by small coupling methods ) 




a: 



DGLAP 



a « 1 



Q 



QCD 



Figure 8. "Phase diagram" of high energy QCD. 
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as we have seen in the previous Section. The partons in the hadronic wave 
function reach the state of saturation. The region of saturation of the struc- 
ture functions is represented in yellow in Fig. ^. The scale Q 2 corresponding 
to the transition to the saturation region is different for different values of 
Bjorken x, increasing with decreasing x (increasing energy s). This scale is 
the saturation scale Q 2 (Y = lnl/x) given by Eq. (|3l|). If Q 2 s (x) > ^qcd 
than the coupling constant inside the saturation region would still be small 
allowing us to perform analytical calculations inside of the saturation regioncj 
by employing Eq. (||). 

It is possible that the recent experimental data obtained at HERA E3 
contains evidence of saturation transition in DIS on a proton at x ~ 10~ 3 — 
10~ 4 and Q 2 f=a 2 — 4GeV 2 , providing an experimental confirmation of the 
theoretical ideas presented above. 

5 Open Questions 

Several question are still open in the field of saturation physics. Here are some 
of the questions related to total cross sections: _ 

1. The large N c limit was instrumental in deriving Eq. (|2^)tiJ. Beyond 
it the systems .of evolution equations in the existing approaches do not seem 
to closc.E3c3n3 Thus the question remains: can one describe the high energy 
asymptotics of QCD beyond the-large-A^ limit? 

2. NLO BFKL correctionstffl seem to be very large andoiegative, pos- 
ing several open questions for the linear evolution equation.Em Would these 
problems be cured if one calculates the NLO correction to the full nonlinear 
evolution? Nonlinear effects become important before the NLO effects and 
may modify the NLO corrections. Work on this project is under way.E3 

3. Finally in this talk we have demonstrated how unitarization happens 
in DIS cross sections via summation of fan diagrams of Fig. |[ This mecha- 
nism would not work for onium-onium scattering cross section, even though 
a large Q~ sho uld exist in the problem. One has to resum the pomeron loop 
diagramsEjEU in order to understand the unitarization of onium-onium cross 
section (and, correspondingly, the proton-proton cross section). Resummation 
of pomeron loops is still an important open question in high energy physics. 
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